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Abstract. The thermodynamics and the dynamics of particle systems with infinite- 
range coupling display several unusual and new features with respect to systems with 
short-range interactions. The Hamiltonian Mean Field (HMF) model represents a 
paradigmatic example of this class of systems. The present study addresses both at- 
tractive and repulsive interactions, with a particular emphasis on the description of 
clustering phenomena from a thermodynamical as well as from a dynamical point of 
view. The observed clustering transition can be first or second order, in the usual 
thermodynamical sense. In the former case, ensemble inequivalence naturally arises 
close to the transition, i.e. canonical and microcanonical ensembles give different re- 
sults. In particular, in the microcanonical ensemble negative specific heat regimes and 
temperature jumps are observed. Moreover, having access to dynamics one can study 
non-equilibrium processes. Among them, the most striking is the emergence of coherent 
structures in the repulsive model, whose formation and dynamics can be studied either 
by using the tools of statistical mechanics or as a manifestation of the solutions of an 
associated Vlasov equation. The chaotic character of the HMF model has been also 
analyzed in terms of its Lyapunov spectrum. 



1 Introduction 

Long-range interactions appear in the domains of gravity jl]|2) and of plasma 
physics Q3| and make the statistical treatment extremely complex. Additional 
features are present in such systems at short distances: the gravitational potential 
is singular at the origin and screening phenomena mask the Coulomb singularity 
in a plasma. This justifies the introduction of simplified toy models that retain 
only the long-range properties of the force, allowing a detailed description of the 
statistical and dynamical behaviors associated to this feature. In this context a 
special role is played by mean-field models, i.e. models where all particles interact 
with the same strength. This constitutes a dramatic reduction of complexity, 
since in such models the spatial coordinates have no role, since each particle 
is equivalent. However, there are several preliminary indications that behaviors 



2 T. Dauxois et al. 



found in mean-field models extend to cases where the two-body potential decays 
at large distances with a power smaller than space dimension Q,|5|,|| ■ 

The Blume-Emery-Griffiths (BEG) mean-field model is discussed in this 
book H and represents an excellent benchmark to discuss relations between 
canonical and microcanonical ensembles. Indeed, this model is exactly solvable 
in both ensembles and is, at the same time, sufficiently rich to display such 
interesting features as negative specific heat and temperature jumps in the mi- 
crocanonical ensemble. Since these effects cannot be present in the canonical en- 
semble, this rigorously proves ensemble inequivalence. However, the BEG model 
has no dynamics and only the thermodynamical behavior can be investigated. 
Moreover, it is a spin model where variables take discrete values. It would there- 
fore be amenable to introduce a model that displays all these interesting ther- 
modynamical effects, but for which one would also dispose of an Hamiltonian 
dynamics with continuous variables, whose equilibrium states could be studied 
both in the canonical and in the microcanonical ensemble. Having access to 
dynamics, one could moreover study non equilibrium features and aspects of 
the microscopic behavior like sensitivity to initial conditions, expressed by the 
Lyapunov spectrum @. Such a model has been introduced in Ref. [|| and has 
been called the Hamiltonian Mean Field (HMF) model. In the simpler version, 
it represents a system of particles moving on a circle, all coupled by an equal 
strength attractive or repulsive cosine interaction. An extension of it to the case 
in which particles move on a 2D torus has been introduced in Ref. and it has 
been quite recently realized that all such models are particular cases of a more 
general Hamiltonian fici|| . 

The HMF model, that we introduce in Section 2, is exactly solvable in the 
canonical ensemble by a Hubbard-Stratonovich transformation. The solution in 
the microcanonical ensemble can be obtained only under certain hypotheses 
that we will discuss in Section 3, but detailed information on the behavior in the 
microcanonical ensemble can be obtained by direct molecular dynamics (MD) 
simulations. The model has first and second order phase transitions and tri- 
critical points. Its rich phase diagram allows to test the presence of ensemble 
inequivalence near canonical first order phase transitions and, indeed, we find 
negative specific heat and temperature jumps in the microcanonical ensemble. 
Having access to dynamics, one can study metastability of out-of-cquilibrium 
states. This is done in Section 4, where we analyze the emergence of a coherent 
structure in the repulsive HMF at low energy. Similar features are also discussed 
in another chapter of this book [[ll| for the attractive case near the second order 
phase transition. Section 5 is devoted to the study of the spectrum of Lyapunov 
exponents. The maximal exponent has a peak near the phase transition P, |l2] , |T3"| 
and vanishes when increasing the number of particles with a universal scaling law 
in the whole high energy disordered phase. In a low energy range the Lyapunov 
spectrum has a thermodynamic limit distribution similar to the one observed 
for systems with short-range interaction 



The Hamiltonian Mean Field Model 



3 



2 The HMF models 



A generic two-body potential in a two dimensional square box of side 2tt with 
periodic boundary conditions, a 2D torus, can be Fourier expanded as 

V(x,y)= jjp exp(ik-T)V(k x ,k y ) . (1) 

A sufficiently rich family of potential functions is obtained if, we restrict to the 
first two momentum shells |A;| = 1 and \k\ = we require that the potential 
is only invariant under discrete rotations by all multiples of 7r/4, and we assume 
that the Fourier coefficients on each shell are the same. This amounts to perform 
a truncation in the Fourier expansion of the potential (Q) , as done in studies of 
spherically symmetric gravitational systems in another chapter of this book |l5| . 
We get 

V(x, y) = a + 6 cos a; + b cosy + c cos x cosy . (2) 

As the constant a is arbitrary and scaling b is equivalent to scale the energy, c 
remains the only free parameter. We consider TV particles interacting through 
the two-body potential V(x,y) and we adopt the Kac prescription [|l6) to scale 
the equal strength coupling among the particles by their number N. This scaling 
allows to perform safely the thermodynamic limit, since both the kinetic and the 
potential energy increase proportional to iVQ. By appropriately redefining the 
constants a — 2e + A, b ~ — e, c = — A in formula (||) and using Kac prescription 
one gets the following potential energy 



+A (1 - cos(a; i - Xj) cos(y 4 - yj)) , (3) 



with (pa, yi) G] — 7r, 7r] x] — tt, tt] representing the coordinates of i-th particle and 
(Px,i,Py,i) its conjugated momentum. The Hamiltonian of the HMF model is 
now the sum of this potential energy with the kinetic energy 




K = \( ; ! * ui 

We get 

Hhmf =K + V a . (5) 

In the following, we will consider model for A = and both e positive 
(attractive case) and negative (repulsive case). The A ^ case will always have 
A > and e = 1. 



1 Kac prescription is, however, unphysical and it would be important to find a viable 
alternative. 
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3 Equilibrium Thermodynamics 

In this section, we discuss the equilibrium thermodynamical results for model (||) 
in the canonical and microcanonical ensembles. Canonical results will be ob- 
tained analytically while, for the microcanonical ones, we will mostly rely on 
molecular dynamics (MD) simulations. 



3.1 Canonical Ensemble for A — 

For pedagogical reasons, we will initially limit our analysis to the case A = 0, 
for which the model reduces to two identical uncoupled systems: one describing 
the evolution of the {xi,p x .i\ variables and the other {yi,p y ,i}- Therefore let us 
rewrite the Hamiltonian associated to one of these two sets of variables, named 
8i in the following. We obtain 

N 2 N 

ffo = £^- + ^£[l- oos(0, - 6j)] =K + V (6) 

i—l — l 

where 6i 6 [— 7r;7r[ and pi are the corresponding momenta. This model can be 
seen as representing particles moving on the unit circle, or as classical JTY"-rotors 
with infinite range couplings. For e > 0, particles attract each other and rotors 
tend to align (ferromagnetic case), while for e < 0, particles repel each other 
and spins tend to anti-align (antiferromagnetic case). At short distances, we can 
either think that particles cross each other or that they collide elastically since 
they have the same mass. 

The physical meaning of this model is even clearer if one introduces the mean 
field vector 

1 N 

M = Me i * = -J2™i (7) 

i=l 

where = (cos 0j, sin 6?,). M and <j> represent the modulus and the phase of the 
order parameter, which specifies the degree of clustering in the particle interpre- 
tation, while it is the magnetization for the XY rotors. Employing this quantity, 
the potential energy can be rewritten as a sum of single particle potentials 

1 N 

V Q = - ^2 v i witn «i = 1 - M cos(8i - 4>) . (8) 

i=l 

It should be noticed that the motion of each particle is coupled to all the others, 
since the mean-field variables M and <f> are determined at each time t by the 
instantaneous positions of all particles. 

The equilibrium results in the canonical ensemble can be obtained from the 
evaluation of the partition function 

Z = [ d N p l d N l exp (-(3H) (9) 
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where (3 — l/(fcsT), with ks the Boltzmann constant and T the temperature. 
The integration domain is extended to the whole phase space. Integrating over 
momenta, one gets: 



Z 



=(ff£ 



d N 9i exp 



-PeN 



(1 -M 2 ) 



(10) 



In order to evaluate this integral, we use the two dimensional Gaussian identity 



exp 



M 2 
2 



OO pOC 



dy exp [-y 2 + • y] 



(11) 



oo J — OO 



where x and y are two-dimensional vectors and fi is positive. We can therefore 
rewrite Eq. (|o|) as 



2 = 'f 



N/2 



exp 



-fieN 



J 



with 



^ />7T />OG />00 

J=- d N 6i / dyexp[-y 2 + V^M-y] 

" J —it J —oo J —oo 



(12) 



(13) 



and /i = (3eN. We use now definition (Q) and exchange the order of the integrals 
in ([jJJ), factorizing the integration over the coordinates of the particles. Intro- 



ducing the rescaled variable y — > yy / 7V/2/3e, one ends up with the following 
expression for J 



J 



N 
2iT(3e 



OO p OO 



dy exp 



oo J —oo 



-iV(^-ln(27T/o(y)) 



(14) 



where I n is the modified Bessel function of order n and y is the modulus of y. 
Finally, integral ( |l4"| ) can be evaluated by employing the saddle point technique 
in the mean-field limit (i.e. for N — > oo). In this limit, the Hclmholtz free energy 
per particle / reads as : 

PS = - ton, ^ = -\ In (|) + f + max - In^))) . (15) 

The maximum condition leads to the consistency equation 



y 



Pe I (y) 



(16) 



For e < 0, there is a unique solution y = 0, which means that the order param- 
eter remains zero and there is no phase transition (see Figs. [l]c, d). The particles 
are all the time homogeneously distributed on the circle and the rotors have zero 
magnetization. On the contrary, in the ferromagnetic case (e > 0), the solution 
y — is unstable for /3 > j3 c = 2. At /3 = /3 C , two stable symmetric solutions ap- 
pear through a pitchfork bifurcation and a discontinuity in the second derivative 
of the free energy is present, indicating a second order phase transition. 
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Fig. 1. Temperature and magnetization as a function of the energy per particle U 
for ^4 = and |e| = 1 in the ferromagnetic (a),(b) and in the antiferromagnetic case 
(c),(d). Symbols refer to MD data for N = 10 2 and 10 3 , while the solid lines refer to the 
canonical prediction obtained analytically in the mean-field limit. The vertical dashed 
line indicates the critical energy in the ferromagnetic model, located at U c = 0.75, 
/3 C = 2. The inset of panel (d) shows the rescaled magnetization M x \/~N . 



These results are confirmed by an analysis of the order parameter^ 

For £ > 0, the magnetization M vanishes continuously at f3 c (see Fig. la,b), 
while it is always identical to zero in the antiferromagnetic case (see Fig. lc,d). 
Since M measures the degree of clustering of the particles, we have for e > a 
transition from a clustered phase when (3 > j3 c to a homogeneous phase when 
f3 < (3 C . We can obtain also the energy per particle 

d(0f) 1 e , „ 

which is reported for |e| = 1 in Fig. [I]. Panels (c) and (d) of Fig. [l| are limited 
to the range U > because in the antiferromagnetic model a non-homogeneous 



2 This is obtained by adding to the Hamiltonian an external field and taking the 
derivative of the free energy with respect to this field, evaluated at zero field. 
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state, a bicluster, can be generated for smaller energies. The emergence of this 
state modifies all thermodynamical and dynamical features as will be discussed 



in section 4.2 



The dynamics of each particle obeys the following pendulum equation of 
motion 

6 i = -Msm(6 i -<t>) , (19) 

where M and <j> have a non trivial time dependence, related to the motion of 
all the other particles in the system. Equation ( |l9| ) has been very successfully 
used to describe several features of particle motion, like for instance trapping 
and untrapping mechanisms There are also numerical indications jg] and 
preliminary theoretical speculations |L7| that taking the mean-field limit before 
the infinite time limit, the time-dependence would disappear and the modulus 
and the phase of the magnetization become constant. This implies, as we will 
discuss in Section 5, that chaotic motion would disappear. The inversion of these 
two limits is also discussed in the contribution by Tsallis et al fll| . 



3.2 Canonical ensemble for A ^ 



As soon as A > 0, the evolution along the two spatial directions is no more de- 
coupled and the system cannot be described in terms of a single order parameter. 
Throughout all this section e = 1. A complete description of the phase diagram 
of the system requires now the introduction of two distinct order parameters: 



M 2 



E 4 sin; 



where Zi 



N ' N 
Xi or yi and; 



M z exp {i<t> z ) 



x±y 



N 



N 



M„ 



= P z exp (iip z 



M and R 



x+y 



- 1 T. — 1 



(20) 

(21) 
- P: 



It can be shown that on average M x 
therefore, we are left with only two or der parameters. 

Following the approach of section 3.1, the canonical equilibrium properties 
can be derived analytically in the mean-field limit M. We obtain 



M 2 + P 2 





In 



G(M, P; A) 



with 



G 



J ds I (M + V2AP cos exp (M cos s) 



(22) 



(23) 



where s is an integration variable. 
The energy per particle reads as 
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Depending on the value of the coupling constant A, the single particle potential 
defined through Va = \ 53j v i changes its shape, inducing the different clustering 
phenomena described below. For small values of A, Vi exhibits a single minimum 
per cell (x,y) € ([— tt, tt], [— 7t,7t]) (see Fig. ||a for A = 1), while for larger values 
of A four minima can coexist in a single cell (see Fig. ||b for A = 4). In the 
former case only one clustered equilibrium state can exist at low temperatures, 
while in the latter case, when all the four minima have the same depth a phase 
with two clusters can emerge, as described in the following. 




Fig. 2. Single particle potential Vi(x, y) for A = l,e = 1 (a) and A = 4, e = 1(b). 



Since the system is ruled by two different order parameters M and P, the 
phase diagram is more complicated than in the A = case and we observe two 
distinct clustered phases. In the very low temperature regime, the system is in 
the clustered phase CPy. the particles have all the same location in a single 
point-like cluster and M w P w 1. In the very large temperature range, the 
system is in a homogeneous phase (HP) with particles uniformly distributed, 
M « P = 0(1/ y/~N). For A > A2 ~ 3.5, an intermediate two-clusters phase 
CP2 appears. In this phase, due to the symmetric location of the two clusters 
in a cell, M ~ 0(l/\/N) while P ~ 0(1) 0. We can gain good insights on the 
transitions by considering the line Tm (resp. Tp) where M (resp. P) vanishes 
and the phase CP\ (resp. CP2) looses its stability (see Fig. || for more details). 

We can therefore identify the following four different scenarios depending on 
the value of A. 

(I) When < A < Ai — 2/5, one observes a continuous transition from the 
phase CPi to HP. The critical line is located at T M = 1/2 (U M = 3/2 + A) and 
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Fig. 3. Canonical phase diagram of model (g) reporting the transition temperatures 
versus the coupling parameter A. The solid (resp. dashed) lines indicate the Tm (resp. 
Tp) lines. The dots the points where the nature of the transitions change. A\, A2 and 
A3 are the threshold coupling constants that determine the transition scenario / — > IV. 



a canonical tricritical point, located at A\ — 2/5, separates the 1st order from 
the second order phase transition regions^ 

(II) When A\ < A < A 2 « 3.5, the transition between CP\ and HP is first 
order with a finite energy jump (latent heat). Inside this range of values of A 
one also finds microcanonical discontinuous transitions (temperature jumps) as 
we will see in Section 3.2. 

(III) When A2 < A < A3 w 5.7, the third phase begins to play a role and 
two successive transitions are observed: first CP\ disappears at Tm via a first 
order transition that gives rise to CP2; then this two-clusters phase gives rise 
to the HP phase via a continuous transition. The critical line associated to this 
transition is T P = A/ '4 (U P = 3A/4 + 1). 

(IV) When A > A3, the transition connecting the two clustered phases, CP\ 
and CP2, becomes second order. 



3 The tricritical point has been identified by finding first the value of P(M, A, (3) that 
minimizes f(M, P; A, (3) and then substituting the solution P in /. This reduces / 
to a function of M only. Now, the standard procedure described also for the BEG 
model JE| can be applied. This consists in finding the value of f3 and A where both 
the M 2 and the M 4 coefficient of the development of the free energy in powers of M 
vanishes 
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3.3 Microcanonical Ensemble 

As we have anticipated, our microcanonical results have been mostly obtained 
via MD simulations, since we cannot estimate easily the microcanonical entropy 
S analytically for the HMF model^j. However, we show below how far we can 
get, starting from the knowledge of the canonical free energy, using Legendre 
transform or inverse Laplace transform techniques. The following derivation is 
indeed valid in general, it does not refer to any specific microscopic model. 

The relation that links the partition function Z(f3, N) to the microcanonical 
phase-space density at energy E — U ■ N 

u(E, N) = J d N Pl d N 9 l 8{E - H) (25) 

is given by 

/•oo 

Z(/3,N) = dE u(E,N) e~ pE , (26) 
Jo 

where the lower limit of the integral (E — 0) corresponds to the energy of the 
ground state of the model. Expression (J2q) can be readily rewritten as 



Z((3,N) = N / dC/exp 



N(-PU+^Hlo(E,N)) 



(27) 



which is evaluated by employing the saddle-point technique in the mean-field 
limit. Employing the definition of entropy per particle in the thermodynamic 
limit 

lnu>(U,N) 



S(U)= lim 

N->oo 



(28) 



_N 

one can obtain the Legendre transform that relates the free energy to the entropy: 

- (3F((3) = max[-/3t/ + S(U)] with (3 = — . (29) 

Since a direct analytical evaluation of the entropy of the HMF model in the 
microcanonical ensemble is not possible, we are rather interested in obtaining 
the entropy from the free energy. This can be done only if the entropy S is a 
concave function of the energy. Then, one can invert (|29j) , getting 

S(U)=mm[(3(U-m)} with U = . (30) 

However, the assumption that S is concave is not true for systems with long range 
interactions near a canonical first order transition, where a "convex intruder" 
of S appears [5|Jl5]|, which gives rise to a negative specific heat regime in the 
microcanonical ensemble^]. For such cases, we have to rely on MD simulations or 
microcanonical Monte-Carlo simulations |19|. 



4 It would be possible using large deviation technique Jig]. 

5 A convex intruder is present also for short range interactions in finite systems, but 
the entropy regains its concave character in the thermodynamic limit 
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An alternative approach to the calculation of the entropy consists in express- 
ing the Dirac S function in Eq. (|25|) by a Laplace transform. One obtains 

u;(E,N) = — d(3e^ E Z(f3) (31) 

^ J-ioo 

where one notices that /3 is imaginary. As the partition function can be estimated 
for our model, we would just have to analytically continue to complex values of (3. 
By performing a rotation to the real axis of the integration contour, one could 
then evaluate lo by saddle point techniques. However, this rotation requires the 
assumption that no singularity is present out of the real axis: this is not in 
general guaranteed. It has been checked numerically for the A — model fis)] , 
allowing to obtain then an explicit expression for u> and confirming ensemble 
equivalence for this case. 

However, when A > At, a first order canonical transition occurs, which im- 
plies the presence of a convex intruder and makes the evaluation of S(U) through 
Eq. [}0| impossible . For these cases, a canonical description is unable to capture 
all the features associated with the phase transition. For the microcanonical 
entropy, we have to rely here on MD simulations. 

The MD simulations have been performed adopting extremely accurate sym- 
plectic integration schemes po|| , with relative energy conservations during the 
runs of order ~ 1CP 6 . It is important to mention that the CPU time required 
by our integration schemes, due to the mean-field nature of the model, increases 
linearly with the number of particles. 

Whenever we observe canonically continuous transitions (i.e. for A < Ai), the 
MD results coincide with those obtained analytically in the canonical ensemble, 
as shown in Fig. [j] for A = 0. The curve T(U) is thus well reproduced from MD 
data, apart from finite iV effects. It has been however observed that starting from 
"water-bags" initial conditions metastable states can occur in the proximity of 
the transition flT|| . 

When A > Ai, discrepancies between the results obtained in the two en- 
sembles are observable in Fig. ||. For example, MD results in the case A = 1, 
reported in Fig. ^(a), differ clearly from canonical ones around the transition, 
exhibiting a regime characterized by a negative specific heat. This feature is 
common to many models with long-range |pl| , p2[ or power-law decaying interac- 
tions Q as well as to finite systems with short-range forces p9fl . However, only 
recently a characterization of all possible microcanonical transitions associated 
to canonically first order ones has been initiated IgJlCj. 

For A slightly above At, the transition is microcanonically continuous, i.e. 
there is no discontinuity in the T—U relation (this regime presumably extends up 
to A ~ 1.2). Before the transition, one observes a negative specific heat regime 
(see Fig. [|a). In addition, as already observed for the Blume-Emery-Griffiths 
model ||, microcanonically discontinuous transitions can be observed in the 
"convex intruder" region. This means that, at the transition energy, temperature 
jumps exist in the thermodynamic limit. A complete physical understanding of 
this phenomenon, which has also been found in gravitational systems p|, has not 
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Fig. 4. Temperature-energy relation in the coexistence region for A=l (a) and A=4 
(b). Lines indicate canonical analytical results, while circles correspond to microcanoni- 
cal MD simulations. Solid thick lines are equilibrium results, solid thin lines metastable 
states and dashed thin lines unstable states. The dash-dotted line is the Maxwell con- 
struction. Figure (a) refers to a first order transition from CP\ to HP, (b) to discontin- 
uous transition connecting the two clustered phases. In (b) the second order transition 
from CP2 to HP associated to the vanishing of P is also shown. The MD results refer 
to model (0) with TV = 5000 averaged over a time t = 10 6 . 




been reached. For A > 1.2, i.e. above the " microcanonical tricritical point", our 
model displays temperature jumps. This situations is shown in Fig. || for A = 2. 

For A > A2, we have again a continuous transition connecting the two clus- 
tered phases CP\ and CP2 . This is the first angular point in the T — U relation 
at U ~ 3.65 in Fig. ||(b). The second angular point at U ~ 4 is the continuous 
transition, connecting CP2 to HP. The transition at lower energy associated to 
the vanishing of M is continuous in the microcanonical ensemble with a negative 
specific heat, while discontinuous in the canonical ensemble; the dash-dotted line 
indicates the transition temperature in the canonical ensemble, derived using the 
Maxwell construction). The second transition associated to the vanishing of P 
is continuous in both ensembles. 



4 Dynamical Properties I: Out-of-equilibrium states 
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Fig. 5. Evolution of the temperature T versus the energy U in the case A = 2. The 
symbols refer to MD results obtained by successively cooling or warming a certain 
initial configuration. Each simulation has been performed at constant total energy and 
refers to a system of N = 4,000 particles integrated for a time t = 10 6 . The solid 
lines are computed in the canonical ensemble and include also unstable and metastable 
cases. The solid vertical line indicating the transition energy has been estimated via a 
Maxwell construction performed in the microcanonical ensemble (for details see |H|). 



4.1 Metastable states 

Around the critical energy, relaxation to equilibrium depends in a very sensitive 
way on the initial conditions adopted. When one starts with out-of-cquilibrium 
initial conditions in the ferromagnetic case, one finds quasi-stationary (i.e. long 
lived) nonequilibrium states. An example is represented by the so-called "water 
bag" initial condition: all the particles are clustered in a single point and the 
momenta are distributed according to a flat distribution of finite width centered 
around zero. These states have a lifetime which increases with the number of 
particles N, and are therefore stationary in the continuum limit. In correspon- 
dence of these metastable states, anomalous diffusion and Levy walks [^3|, long 
living correlations in /i— space |24| and zero Lyapunov exponents |^5[ have been 
found. In addition, these states are far from the equilibrium caloric curve around 
the critical energy, showing a region of negative specific heat and a continuation 
of the high temperature phase (linear T vs U relation) into the low temperature 
one. It is very intriguing that these out-of-equilibrium quasi-stationary states 
indicate a caloric curve very similar to the one found in the region where one 
gets a canonical first order phase transitions, but a continuous microcanonical 
one, as discussed in section |3.3| . In the latter case, however, the corresponding 
states are stationary also at finite N. The coexistence of different states in the 
continuum limit near the critical region is a purely microcanonical effect, and 
arises after the inversion of the t — * oo limit with the N — ► oo one ]24| , p5 11 1. 



14 T. Dauxois et al. 



Similarly, the antiferromagnetic HMF (A — 0, e = —1), where the particles 
interact through repulsive forces presents unexpected dynamical properties in 
the out-of-equilibrium thermodynamics. On the first sight, the thermodynam- 
ics of this model seems to be less interesting since no phase transition occurs 
as discussed above. However, thermodynamical predictions are again in some 
cases in complete disagreement with dynamical results leading in particular to a 
striking localization of energy. This aspect, as we will show below, is of course, 
again, closely related to the long-range character of the interaction and to the 
fact that such a dynamics is chaotic and self-consistent. We mean by this that 
all particles give a contribution to the field acting on each of them. One calls 



this phenomenon, self- consistent chaos 26 . In addition to the toy model that we 



consider here, we do think that similar emergence of structures, but even more 
importantly, similar dynamical stabilization of out-equilibrium states could be 
encountered in other long-range systems, as we briefly describe at the end of this 
section. 



4.2 The dynamical emergence of the bicluster in the 
antiferromagnetic case 

In the antiferromagnetic case, Eq. (^|) with e = — 1, the intriguing properties 
appear in the region of very small energies. To be more specific, if an initial 
state with particles evenly distributed on the circle (i.e. close to the ground state 
predicted by microcanonical or canonical thermodynamics) and with vanishingly 
small momenta is prepared, this initial condition can lead to the formation of 
unstable states. This process, discovered by chance, is now characterized in full 
detail iHHH. 

As shown by Fig. ^, the density of particles is initially homogeneous. However 
a localization of particles do appear at a given time, in two different points, 
symmetrically located with respect to the center of the circle. This localized 
state, that we call bicluster, is however unstable (as shown again by Fig. ^), 
since both clusters are giving rise to two smaller localized groups of particles: 
this is the reason for the appearance of the first chevron. However, also this state 
is unstable, so that the first chevron disappears to give rise again to a localization 
of energy in two points. This state enhances the formation of a chevron with a 
smaller width and this phenomenon repeats until the width of the chevron is so 
small that one does not distinguish anymore its destabilization. Asymptotically 
one gets a density distribution displaying two sharp peaks located at distance tt 
on the circle, a dynamically stable bicluster. 



As we have shown in Ref. 27 , the emergence of the bicluster is the signature 
of shock waves present in the associated hydrodynamical equations. Indeed, we 
found a strikingly good description of the dynamics of the particles by a non 
linear analysis of the associated Vlasov equation, which is mathematically justi- 
fied po|| in the infinite N limit. The physical explanation of this problem can be 
summarized as follows. Once the Hamiltonian has been mapped to the Vlasov 
equation, it is possible to introduce a density p(0,t) and a velocity field v{9,t). 
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Fig. 6. Short-time evolution of the particle density in grey scale: the darker the grey, 
the higher the density. Starting from an initial condition with all the particles evenly 
distributed on the circle, one observes a very rapid concentration of particles, followed 
by the quasi periodic appearance of chevrons, that shrink as time increases. 



Neglecting the dispersion in momentum and relying on usual non linear dynam- 
ics hierarchy of time-scales, one ends up with dynamical equations at different 
orders in a multiscale analysis. The first order corresponds to the linear dynam- 
ics and defines the plasma frequency of order one. However, a second timescale 
appears that is related to the previous one by the relationship r = \/\J t, where 
U is the energy per particle. When one considers initial conditions with a very 
small energy density, the two time scales are very different and clearly distin- 
guishable by considering particle trajectories: a typical trajectory corresponds 
indeed to a very fast motion with a very small amplitude, superimposed to a 
slow motion with a large amplitude. 

This suggests to average over the fast oscillations and leads to the spatially 
forced Burgers equation 

du du 1 . .„„. 

d- T +u ae=-2 sm2e • (32) 

once the average velocity u(8, r) = (v(8, t, tYu is introduced. Due to the absence 



of dissipative or diffusive terms, Equation (32) supports shock waves and these 



can be related to the emergence of the bicluster. By applying the methods of 



characteristics to solve Equation (|32[), one obtains 

d 2 8 1 



, , -sin 20 = . (33) 



dT 2 2 
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which is a pendulum-like equation. 

Fig. shows the trajectories, derived from Eq. (|33|), for particles that are 
initially evenly distributed on the circle. One clearly sees that two shock waves 
appear and lead to an increase of the number of particles around two particular 
sites, which depend on the initial conditions: this two sites correspond to the 
nucleation sites of the bicluster. Because of the absence of a diffusive term, the 
shock wave starts a spiral motion that explains the destabilization of the first 
bicluster and also the existence of the two arms per cluster, i.e. the chevron. 




Fig. 7. Five successive snapshots of the velocity profiles u(8, t) are shown including the 
initial state when all particles are uniformly distributed in space with a small velocity 
dispersion. 



This dynamical analysis allows an even more precise description, since the 
methods of characteristics show that the trajectories correspond to the motion 
of particles in the double well periodic potential V(6,t) = — (cos2#)/4. This 
potential is of mean field origin, since it represents the effect of all interacting 
particles. Therefore, the particles will have an oscillatory motion in one of the 
two wells. One understands thus that particles starting close to the minimum 
will collapse at the same time, whereas a particle starting farther will have a 
larger oscillation period. This is what is shown in Fig. || where trajectories are 
presented for different starting positions. One sees that the period of recurrence 
of the chevrons corresponds to half the period of oscillations of the particles 
close to the minimum; this fact is a direct consequence of the isochronism of the 
approximate harmonic potential close to the minimum. The description could 
even go one step further by computing the caustics, corresponding to the en- 
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velops of the characteristics. They are shown in Fig. ^ and testify the striking 
agreement between this description and the real trajectories: the chevrons of 
Fig. ^ correspond to the caustics. 




Fig. 8. Superposition of the caustics (thick full lines) over the characteristics (dotted 
lines) of particles that are initially evenly distributed between —it/2 and n/2. 



4.3 Thermodynamical predictions versus dynamical stabilization 

If the above description is shown to be particularly accurate, it does not explain 
why this state is thermodynamically preferred over others. Indeed, as shown be- 
fore, thermodynamics predicts that the only equilibrium state is homogeneous. 
This result has been discussed in Section 3.1, where we have proved that magne- 
tization is zero at all energies for the antifcrromagnetic model and this has been 
also confirmed by Monte Carlo simulations |2q| . However, since MD simulations 
are performed at constant energy, it is important to derive analytically the most 
probable state in the microcanonical ensemble. Since this model does not present 
ensemble inequivalence, we can obtain the microcanonical results by employing 
the inverse Laplace transform ( |3l| ) of the canonical partition function. The mi- 
crocanonical solution confirms that the maximal entropy state is homogeneous 
on the circle. It is therefore essential to see why the bicluster state, predicted to 
be thermodynamically unstable is instead dynamically stable. 

The underlying reason rely on the existence of the two very different timesca- 
les and the idea is again to average over the very fast one. Instead of using 
the classical asymptotic expansion on the equation of motions, it is much more 
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appropriate to develop an adiabatic approximation which leads to an effective 
Hamiltonian that describes very well the long time dynamics. Doing statistical 
mechanics of this averaged problem, one predicts the presence of the bicluter. 

The theory that we have developed relies on an application of adiabatic 
theory, which in the case of the HMF model is rather elaborate and needs lengthy 
calculations j2?|] that we will not present here. An alternative, but less powerful, 
method to derive similar results has been used in Ref. p9fl . On the contrary, we 
would like to present a qualitative explanations of this phenomenon, using a nice 
(but even too simple !) analogy. 

This stabilization of unstable states can be described using the analogy with 
the inverted pendulum, where the vertical unstable equilibrium position can be 
made stable by the application of a small oscillating force. One considers a rigid 
rod free to rotate in a vertical plane and whose point-of-support is vibrated ver- 
tically as shown by Fig. ||a. If the support oscillates vertically above a certain 
frequency, one discovers the remarkable property that the vertical position with 
the center of mass above its support point is stable (Fig. ^b) . This problem, dis- 
cussed initially by Kapitza |3l| , has strong similarities with the present problem 
and allows a very simplified presentation of the averaging technique we have 
used. 




a cos(cjt) 



a cos(cjt) 



Fig. 9. Schematic picture of the inverted pendulum. 



The equation of motion of the vibrating pendulum is 
d 2 I ' ui 2 \ 

— 2+ i ^§ -acosi) sin0 = O , (34) 

where loq is the proper linear frequency of the pendulum, lu the driving frequency 
of the support and a the amplitude of excitation. Introducing a small parameter 
e = luq/lo, one sees that ( |34| ) derives from the Lagrangian 

1 / dO \ 2 

£ = -f— J + (e 2 - acost) cos6> . (35) 

Here the two frequencies u> and wo define two different time scales, in close 
analogy with the HMF model. Using the small parameter to renormalize the 
amplitude of the excitation as A — a/e, and choosing the ansatz 9 = 8q(t) + 
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e9\(t, e) where r = et, the Lagrangian equations for the function 6\, leads to the 
solution B\ = — A cost sin#o(T). This result not only simplifies the above ansatz, 
but more importantly suggest to average the Lagrangian on the fast variable t 
to obtain an effective Lagrangian £ e // = {C)t = \ (^57) + Veff , where the 
averaged potential is found to be 

A 2 

V eff = - cos 6»o - — cos 26> + Cste . (36) 
8 

It is now straightforward to show that the inverted position would be stable if 
A 2 > 2, i.e. if auj > \2ujq. As the excitation amplitude a is usually small, this 
condition emphasizes that the two time scales should be clearly different, for the 
inverted position to be stable. 

The procedure for the HMF model is analogous, but of course it implies a 
series of tedious calculations. Since we would like to limit here to a pedagogical 
presentation, we will skip such details that can be found in Ref. |^7| . It is however 
important to emphasize that the potential energy in the HMF model is self- 
consistently determined and depend on the position of all particles. The magic 
and the beauty is that, even if this is the potential energy of N particles, it is 
possible to compute the statistical mechanics of the new effective Hamiltonian, 
derived directly from the effective Lagrangian via the Legendre transform. The 
main result is that the out-of-equilibrium state, (i.e. the bicluster shown in Fig. ^) 
corresponds to a statistical equilibrium of the effective mean-field dynamics. No 
external drive is present in this for the inverted pendulum, but the time 

dependence of the mean field plays the role of the external drive. 

The HMF model represents presumably the simplest A-body system where 
out-of-equilibrium dynamically stabilized states can be observed and explained 
in detail. However, we believe that several systems with long range interactions 
should exhibit behaviours similar to the ones we have observed here. Moreover, 
this model represents a paradigmatic example for other systems exhibiting non- 
linear interactions of rapid oscillations and a slower global motion. One of this 
is the piston problem ]32]| : averaging techniques have been applied to the fast 
motion of gas particles in a piston which itself has a slow motion f33|| . Examples 
can also be found in applied physics as for instance wave-particles interaction 
in plasma physics |3~if| , or the interaction of fast inertia gravity waves with the 
vortical motion for the rotating Shallow Water model |35fl . 

5 Dynamical Properties II: Lyapunov exponents 

In this section, we discuss the chaotic features of the microscopic dynamics of the 
HMF model. We mainly concentrate on the A — case, presenting in details the 
behavior of the Lyapunov exponents and the Kolmogorov-Sinai entropy both for 
ferromagnetic and antiferromagnetic interactions. We also briefly discuss some 
peculiar mechanisms of chaos in the A ^ case. The original motivation for 
the study of the chaotic properties of the HMF was to investigate the relation 
between phase transitions, which are macroscopic phenomena, and microscopic 
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dynamics (see Ref. |56| for a review) with the purpose of finding dynamical 
signatures of phase transitions |57j. Moreover, we wanted to check the scaling 
properties with the number of particles of the Lyapunov spectrum [ fl4"| in the 
presence of long-range interactions. 

5.1 The A — case 

In the A — case (see formulae (p|)-(|l9|)), the Hamiltonian equations of motion 
are 

0i = Pi (37) 
Pi = -Msin(0i-(£) . (38) 

The Largest Lyapunov Exponent (LLE) is defined as the limit 

A 1= lim -ln^rr (39) 
t^oo t d(0) K ' 



where d(t) = y^^fffli) 2 + (Spi) 2 is the Euclidean norm of the infinitesimal 
disturbance at time t. Therefore, in order to obtain the time evolution of d(t), 
one has to integrate also the linearized equations of motion along the reference 
orbit 

i sei=5pi ' j t Spi =-£flS** ' (40) 

3 J 

where the diagonal and off-diagonal terms of the Hessian are 

=Mcos(0 i -*)-l (41) 
d 2 V 1 



d6id6j N 



coaiOi-Oj) , i+3 ■ (42) 



To calculate the largest Lyapunov exponent we have used the standard method 
by Benettin et al |3q ]. In Fig. (^0|), we report the results obtained for four 
different sizes of the system (ranging from N — 100 to N = 20, 000). 

In panel (a), we plot the largest Lyapunov exponent as a function of U. As 
expected, Ai vanishes in the limit of very small and very large energies, where the 
system is quasi-integrable. Indeed, the Hamiltonian reduces to weakly coupled 
harmonic oscillators in the former case or to free rotators in the latter. For 
U < 0.2, Ai is small and has no iV-dependence. Then it changes abruptly and a 
region of "strong chaos" begins. It was observed || that between U = 0.2 and 
U = 0.3, a different dynamical regime sets in and particles start to evaporate 
from the main cluster, in analogy with what was reported in other models [3^j39|| . 
In the region of strong chaoticity, we observe a pronounced peak already for 
N = 100 @. The peak persists and becomes broader for N = 20,000. The 
location of the peak is slightly below the critical energy and depends weakly 
on N. 
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Fig. 10. Largest Lyapunov exponent LLE and kinetic energy fluctuations E — ok /v~N 
as a function of U in the A — ferromagnetic case for different N sizes. The theoretical 
curve is shown as a full line, see text. 



In panel (b), we report the standard deviation of the kinetic energy per 
particle £ computed from 



OK_ = y/(K*) - (K) 2 



N 



(43) 



where (•) indicates the time average. The theoretical prediction for E, which is 
also reported in Fig.[l0[ is |4l|J~ 



T 

71 1 



1 - 2M 



V dT 



(44) 



where M(T) is computed in the canonical ensemble. Finite size effects are also 
present for the kinetic energy fluctuations, especially for U > U c , but in general 
there is a good agreement with the theoretical formula, although the experi- 
mental points in Fig. flQb lies systematically below it. The figure emphasizes 
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that the behavior of the Lyapunov exponent is strikingly correlated with S: 
in correspondence to the peak in the LLE, we observe also a sharp maximum 
of the kinetic energy fluctuations. The relation between the chaotic properties 
and the thermodynamics of the system, namely the presence of a critical point, 
can be made even more quantitative. An analytical formula, relating (in the 
A = model) the LLE to the second order phase transition undergone by the 
system, has been obtained Q by means of the geometrical approach developed 
in Refs. BJ3^]. Using a reformulation of Hamiltonian dynamics in the language 
of Ricmannian geometry, they have found a general analytical expression for the 
LLE of a Hamiltonian many-body system in terms of two quantities: the average 
Qq and the variance an of the Ricci curvature = AV/N = 1/N 
where V is the potential energy and qi are the coordinates of the system. Since 
in the particular case of the HMF model, we have 

t=i 1 

the two quantities SIq and an can be expressed in terms of average values and 
fluctuations either of M 2 or of the kinetic energy K 

« =(M 2 )-i = l W + (l-2 f /)-i 

a 2 n = Na 2 M2 = lo£ . (46) 



The formula obtained by Firpo 41 relates the LLE, a characteristic dynamical 
quantity, to thermodynamical quantities like (M 2 ) and a M 2, or (K) and ax, 
which characterize the macroscopic phase transition. For moderately small values 
of U , an approximation of the formula gives 

A oc -^L = S . (47) 

Vn 

This is in agreement with the proportionality between LLE and fluctuations of 
the kinetic energy found numerically in Fig. ([To|) . This implies also a connection 
between the LLE and the specific heat, another quantity which is directly related 
to the kinetic energy fluctuation. In fact the specific heat can be obtained from 
S by means of the Lebowitz-Percus-Verlet formula E3 



2I T 



(48) 



In Fig. [ll], we report the numerical results for the specific heat as a function 
of U for a system made of N = 500 particles, and we compare them with the 
theoretical estimate. 

In the HMF model and for a rather moderate size of the systems, it is possible 
to calculate not only the LLE but all the Lyapunov exponents, and from them the 
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Fig. 11. Specific heat as a function of U in the ^4 = ferromagnetic case. The numerical 
simulation at equilibrium for a system with N = 500 is compared with the expected 
theoretical result (|l8|), 

Kolmogorov-Sinai entropy. We give first a succinct definition of the spectrum of 
Lyapunov exponents (for more details see JtJ ) . Once the 2N-dimensional tangent 
vector z = {§6\, • • • , 86 n, Sp±, ■ ■ ■ ,Spn) is defined, with its dynamics given by 
Eqs. (pOj), one can formally integrate the motion in tangent space up to time t, 
since the equations are linear, 

z(t) = J*z(0) , (49) 

where J* is a 27V x 2N matrix that depends on time through the orbit 0i(t),pi(t). 
The first k exponents of the spectrum Ai, . . . , A&, which are ordered from the 
maximal to the minimal, are then given by 

(Ai+... + A fe )= lim llnTrJ*(4)* , (50) 

t — >oo Zt 

where J l k is the matrix ((J^)* its transpose) induced by J* that acts on the 
exterior product of k vectors in the tangent space z± A . . . A z&. The spectrum 
extends up to k = 2N and in our Hamiltonian system obeys the pairing rule 

Xi = -X 2 N+i-i for 1 < i < 2N . (51) 

The numerical evaluation of the spectrum of the Lyapunov exponents is 
a heavy computational task, in particular for the necessity to perform Gram- 
Schmidt orthonormalizations of the Lyapunov eigenvectors in order to maintain 
them mutually orthogonal during the time evolution. We have been able to 
compute the complete Lyapunov spectrum for system sizes up to N = 100. In 
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Fig. |l2L we report the positive part of the spectrum for different system sizes 
and an energy U — 0.1 inside the weakly chaotic region. The negative part of 
the spectrum is symmetric due to the pairing rule (|5l"|). The limit distribution 
A(x), suggested for short range interactions, 

X(x) = lira X xN (N) , (52) 

that is obtained by plotting Ai vs. i/N and letting N going to infinity, is found 
also here for the ./V values that we have been able to explore. At higher energies, 
this scaling is not valid and a size-dependence is present G3] . 
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Fig. 12. Scaling of the positive part of the spectrum of Lyapunov exponents in the 
A=0 ferromagnetic case for U = 0.1. 

The Kolmogorov-Sinai (K-S) entropy is, according to Pesin's formula 0, the 
sum of the positive Lyapunov exponents. In Fig. |l3|, we plot the entropy density 
Sks/N as a function of U for different systems sizes. 

As for the LLE, Sks/N shows a peak near the critical energy, a fast con- 
vergence to a limiting value as N increases in the small energy limit, and a 
slow convergence to zero for U > U c . A comparison of the ferromagnetic and 
antiferromagnetic cases is reported in Fig. [l4|. Here, for N=100, we plot as a 
function of U the LLE and the Kolmogorov-Sinai entropy per particle Sks/N. 
In both the ferromagnetic and antiferromagnetic cases, the system is integrable 
in the limits of small and large energies. The main difference between the fer- 
romagnetic and the antiferromagnetic model appears at intermediate energies. 
In fact, although both cases are chaotic (LLE and Sks/N are positive), in the 
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Fig. 13. Sks/N as a function of U in the A = ferromagnetic case. Numerical calcu- 
lations for different systems sizes ranging from TV = 10 to TV = 100 are shown. The 
dashed line indicates the critical energy. 



ferromagnetic system one observes a well defined peak just below the critical 
energy, because the dynamics feels the presence of the phase transition. On the 
other hand, a smoother curve is observed in the antiferromagnetic case. In the 
low energy regime, it is possible to work out jL2fl a simple estimate Ai oc VU, 
which is fully confirmed for the ferromagnetic case in Fig. |l5| for different sys- 
tem sizes in the range TV = 100, 20000. The same scaling law is also valid in the 
antiferromagnetic case and for the A = 1 Q HMF model. 

At variance with the TV-independent behavior observed at small energy, 
strong finite size effects are present above the critical energy in the ferromagnetic 
case and for all energies for the antiferromagnetic case. In Fig. [L6|(a), we show 
that the LLE is positive and TV-independent below the transition (see the values 
U = 0.4,0.5), while it goes to zero with TV above. We also report in the same 
figure a calculation of the LLE using a random distribution of particle positions 
6i on the circle in Eq. ( ^fj| ) for the tangent vector. The agreement between the de- 
terministic estimate and this random matrix calculation is very good. The LLE 
scales as TV - 3 , as indicated by the fit reported in the figure. This agreement can 
be explained by means of an analytical result obtained for the LLE of product 
of random matrices |fl6|| . If the elements of the symplectic random matrix have 
zero mean, the LLE scales with the power 2/3 of the perturbation. In our case, 
the latter condition is satisfied and the perturbation is the magnetization M. 
Since M scales as TV" 2 ; we get the right scaling of Ai with TV. This proves that 
the system is integrable for U > U c as TV — > 00. This result is also confirmed by 
the analytical calculations of Ref. |Q and, more recently, of Ref. [ ^p[ . 
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Fig. 14. LLE and Sks/N as a function of U in the A = ferromagnetic and antifer- 
romagnetic cases for N=100. 



In the antiferromagnetic case, the LLE goes to zero with system size as N~i 
for all values of U . 

Interesting scaling laws have also been found for the Kolmogorov- Sinai en- 
tropy in the ferromagnetic case: at small energies Sks/N oc U 3/4 with no size 
dependence, and Sks/N oc N^ 1 / 5 for overcritical energy densities. The latter 
behavior has been found also in other models [^Tj. Concluding this section we 
would like to stress that the finite value of chaotic measures close to the critical 
point is strongly related to kinetic energy fluctuations and can be considered as 
a microscopic dynamical indication of the macroscopic equilibrium phase transi- 
tion. This connection has been found also in other models and seems to be quite 
general @0g|,|§ 

The behavior of the HMF model as a function of the range of the interac- 
tion ]l3|,p|J^| and the dynamical features before equilibration |^,^5) is discussed 
in a separate chapter of this volume in connection to Tsallis nonextensive ther- 
modynamics |Q. 



5.2 Mechanisms of chaos in the 4^0 case 

For A = 0, the origin of chaos is related to the non time-dependence of Eq. ([To|) , 
since it is obvious that if the phase <p and the magnetization M would become 
constant the dynamics of the system will reduce to that of an integrable system. 
There are indeed preliminary indications fl7| that in the mean-field limit N — > 
oo, M and 4> will become constant and A — ► . It should be noticed that this 
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is true if the mean-field limit is taken before the limit t — > oo in the definition 
of the maximal Lyapunov exponent, and numerical indications were reported 
in Ref. When A > we expect a quite different situation: indeed, even 
assuming that in the mean-field limit M and P and their respective phases will 
become constant, the dynamics will eventually take place in a 4-dimensional 
phase space and chaos can in principle be observed. 

As already shown in [||, for A = 1 and e = 1, two different mechanisms of 
chaos are present in the system for U < U c : one acting on the particles trapped 
in the potential and another one felt by the particles moving in proximity of the 
separatrix. This second mechanism is well known and is related to the presence of 
a chaotic layer situated around the separatrix. The origin of the first mechanism 
is less clear, but presumably related to the erratic motion of the minimum of 
the potential well, i.e. to the time-dependent character of the equations ruling 
the dynamics of the single particle. Indications in this direction can be found by 
performing the following numerical experiment. Let us prepare a system with 
N = 200 and U = 0.87 (the critical energy is in this case U c ~ 2) with a 
Maxwellian velocity distribution and with all particles in a single cluster. 

For an integration time t < 2 x 10 6 , the Lyapunov exponent has a value 
A ~ 0.13. But when at time t ~ 2 x 10 6 , one particle escapes from the cluster, 
its value almost doubles (see Fig. [l7]). The escaping of the particle from the 
cluster is associated to a decrease of the magnetization M and of the kinetic 
energy K. This last effect is related to the negative specific heat regime: the 
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Fig. 15. Scaling properties of the LLE at low energies in the A = ferromagnetic case. 
No N-dependence is observed for U < 0.2. The dashed line indicates a power-law U 1 ^ 2 . 
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Fig. 16. Scaling of the LLE vs TV for the A = ferromagnetic and antiferromagnetic 
cases at various energies, see text. A power-law decreasing as iV -1 / 3 of the LLE is 
observed for overcritical energies in the ferromagnetic case and for all energies in the 
antiferromagnetic one. See text for further details. 



potential energy Va is minimal when all the particles are trapped, if one escapes 
then Va increases and due to the energy conservation K decreases. As a matter 
of fact, we can identify a "strong" chaos felt from the particles approaching the 
separatrix and a "weak" chaos associated to the orbits trapped in the potential 
well. We believe that the latter mechanism of chaotization should disappear (in 
analogy with the ^4 = case) when the mean-field limit is taken before the 
t — * oo limit. Therefore we expect that for N — > oo the only source of chaotic 
behaviour should be related to the chaotic sea located around the separatrix. As 
already noticed in Ref. || , the degree of chaotization of a given system depends 
strongly on the initial condition (in particular in the mean- field limit). In the 
latter limit, for initial condition prepared in a clustered configuration, we expect 
that A = 0, until one particle will escape from the cluster. 

6 Conclusions 

We have discussed the dynamical properties of the Hamiltonian Mean Field 
model in connection to ita thermodynamics. This apparently simple class of 
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Fig. 17. Time evolution of the Largest Lyapunov Exponent (LLE) Ai, of the magne- 
tization M, and of the kinetic energy K are shown for a clustered initial condition for 
the model with A = l and s = 1 at U = 0.87 and with N = 200. 



models has revealed a very rich and interesting variety of behaviours. Inequiv- 
alence of ensembles, negative specific heat, metastable dynamical states and 
chaotic dynamics are only some among them. During the past years these mod- 
els have been of great help in understanding the connection between dynamics 
and thermodynamics when long-range interactions are present. Such kind of in- 
vestigation is of extreme importance for self-gravitating systems and plasmas, 
but also for phase transitions in finite systems, such as atomic clusters or nuclei, 
and for the foundation of statistical mechanics. Several progresses have been 
done during these years. This contribution, although not exhaustive, is an ef- 
fort to summarize some of the main results achieved so far. We believe that the 
problems which are still not understood will be hopefully clarified in the near 
future within a general theoretical framework. We list three important open 
questions that we believe can be reasonably addresses: the exact solution of the 
model in the microcanonical ensemble; the full characterization of the out-of- 
equilibrium states close to phase transitions; the clarification of the scaling laws 
of the maximal Lyapunov exponent and of the Lyapunov spectrum; the study 
of the single-particle diffusive motion p3f in the various non-equilibrium and 
equilibrium regimes. 
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